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Abstract :Using a previously derived QCD effective hamiltonian we find the masses 
of heavy quarkonia states. Non perturbative effects are included through temperature 
dependent gluonic condensates. We find that even a moderate change in these conden- 
sates in a hot hadronic environment (below the deconfining transition) is sufficient to 
significantly change the heavy meson masses. 



The study of hadronic matter at high temperatures and densities has direct rele- 
vance for heavy-ion experiments. Apart from the search of quark gluon plasma, serious 
attention has been given to signatures of a hot system composed by hadrons below the 
deconfming transition. Among these signatures, changes in the masses due to medium 
effects play a major role. They have been extensively studied with the Nambu-Jona- 
Lasinio model [1], with non-relativistic potential models [2] , with QCD sum- rules [3] 
and in lattice QCD [4]. 

The experimental detection of medium effects on the masses is a very difficult 
problem. On the other hand, from the theoretical point of view the situation is not 
quite clear. Some calculations predict a significant decrease (~ 400 MeV) of the rj c , 
J/ip and tp' masses. Some others suggest that the masses stay constant. We will argue 
that they may increase. 

The approach to this problem adopted in the present work is in many aspects 
similar to QCD sum-rules at finite temperature. In particular, our results for the heavy 
quarkonium spectrum depend primarily on perturbation theory (which is related to the 
value of a s ) and on the gluon condensate. The same conclusion is found in ref. [3]. The 
main difference between this work and the above mentioned spectrum calculations is 
that we give special emphasys to vacuum changes with temperature. 

It has been known since the late seventies that the QCD (physical) vacuum is full of 
soft gluons or, equivalently, contains chromoelectric (E) and magnetic (B) fields. Such 
a state, \ fl), has lower energy than a state without any fields, the perturbative vacuum, 
|0). This picture is supported by the existence of non- vanishing gluon condensates, i.e., 

(n|^^F"^) = ^Vo (i) 

7T 

where is the usual QCD field tensor. 
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At increasing temperatures, lattice calculations predict a phase transition to a de- 
confined phase. In terms of the vacuum state this transition can be interpreted as 
the passage from the non-perturbative (or physical) to the perturbative vacuum, i.e., 
\Q) — > |0). Roughly speaking, one can say that the background soft gluon fields would 
then "boil away", implying that <fi 2 — > 0. In fact some specific lattice calculations [5] 
suggest that even above the deconfining phase transition there may be non-vanishing 
gluon condensates, but they are not yet conclusive and therefore from lattice simula- 
tions we cannot extract the temperature dependence of the gluon condensate over a 
wide range of temperatures. Model calculations using chiral perturbation theory [6] or 
dilaton fields [7] come to the conclusion that <fi 2 stays almost constant until the phase 
transition temperature T c and then starts to drop faster going to zero very slowly. 

We will now shortly describe the basic elements of our "effective QCD" applied 
to the study of heavy quarkonium. Our meson states incorporate the background 
soft gluon fields . Together with the usual quark anti-quark states, \qq), we will 
have also \qEq) or \qB~q). In the first case the quark and anti-quark are in a color 
singlet representation and we call it a singlet state. In the second case, because of the 
interaction with the vacuum, the combination qEq (or qBq . . .) is a color singlet but 
the quark anti-quark pair is in a color octect representation. We call therefore these 
states octect states. 

In order to ensure the gauge invariance of our calculations we will make our meson 
states gauge invariant by construction. This can be done with the help of the color 
transport operator 



The path ordered exponential (denoted by P exp) of the background gluon field A^(x) 
transports along a straight line the color index b (of the fundamental representation of 
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the gauge group) at position x\ to index a at position x 2 - This operator was introduced 
in this context by Schaden and Glazek [8] and extensively used later on by Nunes [9] . It 
is easy to show that canonical anti-commutation relations for the quark and anti-quark 
operators imply that the singlet meson basis states defined below are orthonormal. 

Considering all that was said above we can write now our heavy meson basis states 
explicitly. For simplicity we restrict ourselves to the pseudoscalar mesons, which are 
then 

l*> = 4e*4 q (^) y;„,(^..r;) ,•;;(/,) < (3) 

v 6 ab . a 



\B) = £ -^uf^^-S^Tfc^fOt^^Oin) (4) 

a,o,c;a,p ^ 



l^i) = E ^«i Q (^)^-(x 2 -fi)r 6c (f 2 ,f 1 )^(f 1 )|n> (5) 
\e 2 ) = ]T i ^«f , (^)^-(^x(f 2 -f 1 ))r fc (f 2 ,f 1 )«f(f 1 )|n> (6) 

a,b,c; a ,p V27T0 

The pseudoscalar meson rj can be well represented by a linear combination of the 
above basis states 

\V)= E / lM2,l)|2,l) M . (7) 

M=S,E 1 ,E 2 ,B Jl ' 2 

The interaction between heavy quarks is described by the QCD Lagrangian with 
two simplifying approximations: a) non relativistic limit with the inverse heavy quark 
mass expansion up to first order and b) separation of the gluon fields into classical 
background nonperturbative fields (which will later give rise to the condensates) and 
quantum high momentum fields. Expansions involving both types of fields will include 
only lower order terms because we will consider only the lowest order gluon condensates 
and also because higher powers of the quantum high momentum fields will couple only 
in the perturbative regime and can be neglected since a s is small for high momentum 
couplings. 



Denoting the quark fields by tp and gluon fields by the QCD Lagrangian is 
written as 

C QCD = -\F^Fr + ^(t^ + gT a %)^-m^ . (8) 
We make then a Foldy-Wouthuysen transformation in the quark fields 



ip — > exp (27 • D/2rn) ip 
ip — > ip exp (—ij-D/ 2rnJ 



(9) 



obtaining a non-relativistic Lagrangian 

Cnrqcd = ~ \ F; u Ff + i> (V A) -™)^ + ^D 2 i) + ±-i>gt-Bi) (fO) 

where D = d — igV and E = ^ ^ ^ ^ does not couple upper and lower spinor 
components. We next separate the gluon field in background (A^) and quantum 
fields 

v* = a; + q; . 

We choose the Coulomb background gauge for the quantum fields 

AQ l = (12) 

where D^Q U = d^Q va + gf abc A^ b Q„ c . The background fields are defined in a modified 
Schwinger gauge [10] 



A) = -\F)^ ■ A\ = -F b m x* . (13) 



The field A^ is treated as an external field and therefore satisfies the equation of 
motion D^F^ = 0, where F^ v is the background field strength which is assumed to be 
practically constant over the extent of the heavy meson. 

We next expand the nonrelativistic Lagrangian only to second order in the quan- 
tum fields and subsequently integrate them out in favour of an effective (coulombic) 
interaction. These calculations have been carried out in more detail in ref. [9,11] and 



will not be presented here. It is important to mention that during the calculations 
retardation effects have been neglected (this instantaneous approximation should be 
correct up to order 1/m) and that matrix elements of have been parametrized, 
producing terms proportional to 2 . The resulting effective Hamiltonian was presented 
in ref. [9,11] and is still complicated. We have numerically diagonalized it in the basis 
(3 - 6) and found the solution of the resulting set of coupled differential equations for 
the wave functions. We have then checked that for the description of the low lying 
states of the spectrum it is enough to keep the terms of order (l/m)° plus the kinetic 
energy terms. The effective Hamiltonian can be finally written as 

H = J d 3 x ^u\x) mu{x) + v (x) mv\x) — - F^ u 

— u\x ) T A g Ef Xi u(x ) — v(x ) T A g Ef Xi v\x ) 

V 2 V 2 

- uHx ) — -u(x ) - v(x ) — - vHx ) 

K J 2m K J K J 2m K ' 



f 1 A 

+ a s J d 3 yu\x)T A u(x) - v(y)T v\y) 



(14) 



where the second line corresponds to the "Stark effect" discussed by Leutwyler [14] and 
Voloshin [13]. Here u(x) and v(x) denote the annihilation operators for a quark and 
antiquark of mass m respectively whose spin and color indices have been suppressed, 
r = x — y and T A , T A are the Hermitian generators of the SU (3) color Lie-algebra in 
the 3 and 3 representations respectively. 
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Diagonalysing our effective Hamiltonian in the basis (3-6) only the states S and Ei 
couple and we obtain the following set of coupled differential equations for the wave 
functions. 

(15a) 



n „ 1 d 2 4 a s 
2m -E + C--— - - ^ 

m ar A 6 r 



S(r) = -^E l{ r) 



2m-E--[ — 
m \ or z 



4 +^ 



3^2 



(156) 



r* I 6 

where is the mass eigenvalue of the quarkonium and m is the mass of the constituent 
quarks. The functions S(r) and Ei(r) are related to the wave function components 
in the expansion (7) via 

1 E 

*l>s(r) = -S(r) ; ip El (r) = -± (r) . 

In order to include the scale dependence (or distance dependence) of a s we use in 
eq. (15) 

a =a(r)= + ^ + 53/75 - 

S s[ ) b f(rf + f(r) 625 f(r) J 



where 



m=ln[ 7r r 5 + b] 



1 

6 = 25/3 



This is the result of the two loop calculation given in ref. [12]. In the above expression 
7s = 0.5772, A = 200 MeV and 6 = 20. Among the matrix elements we find some 
involving the energy of the background fields. In particular we find 

(S\ J d 3 x (-If^F"") \S) = (n\F 2 \n) = C = C <P 2 (16a) 

(E ± \ J d 3 x (-^F^F^ \E^ ~ (Q\EF 2 E\Q) = (166) 

where (ft 2 is the gluon condensate defined in (1) , Co is a positive constant and C is the 
energy appearing in eq. (15a). The last line follows from the assumption that higher 
order vacuum expectation values of the background fields are zero. 



As it can be seen there is a splitting between "singlet" and "octet" states given by C. 
Since 2 is negative C will be also negative and the singlet states have a negative energy 
with respect to the octet states. The constant factor Co was fixed by reproducing the 
observed energy levels of the groundstate and first excited charmonium and bottomium 
states at zero temperature. In the (l/m)° approximation pseudoscalar and vector 
mesons are degenerate and our calculations are valid for the J /ip , ip' ,T and T' states. 
For an energy splitting between ip' and J/ip of 600 MeV we obtain C = —737 MeV. 

In order to investigate the temperature dependence of our results we consider the 
temperature dependence of the gluon condensates. In view of the existing estimates of 
this dependence we parametrize it in the folowing way : 



where 0q is the value of the condensate at T = ( 0q = 360MeV 4 ) and T is some 
critical temperature at which 2 = 0. T might be much larger than the deconfining 
transition temperature. With T = 250 MeV eq. (17) interpolates the results compiled 
in ref. [6]. Inserting eq. (17) into (15) and solving it for several values of T between 
and T we find the wave functions and masses of the fundamental and first excited 
states at different temperatures. The results for the masses are shown in figure 1. 
The quark masses were taken to be m c = 1640 MeV and rrib = 4800 MeV. The first 
interesting aspect in figure 1 is that the masses are increasing with temperature. We 
understand this behaviour in the following way : any physical state considered is a 
mixture of singlet and octet components but for the low- lying states , such as J/ip , 
ifj', T and T' the singlet component is largely dominant. In the first of equations 15 
(15a) , the constant C has the effect of shifting the energy of the state to smaller values 
(since it is negative). The temperature dependence of the condensate implies that 




(17) 



C = C{T) = C o 2 (T) 
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As 4> 2 {T) decreases with temperature , so does C and the energy levels of the states are 
shifted to larger values. Roughly speaking we " boil away " the physical vacuum and 
raise the energy of the states. A similar effect occurs in the context of thermodynamics 
of quarks and hadrons during the deconfining transition, where a suppression of the 
physical vacuum brings an additional term (in the simple bag model language this is 
the bag constant) to the quark energy density. 

Our conclusion about the behaviour of heavy quarkonium masses with temperature 
is in contradiction with those of ref. [1] and [2] but in agreement with the lattice 
simulations performed in ref. [4]. 

The second interesting aspect in fig. 1 is the sudden disappearence of the ift' line 
much below the critical temperature. Excited states (specially light ones) are less 
tightly bounded and therefore sensitive to subtle changes in the confining potential. 
Moreover they contain larger octet components , which vary rapidly with the con- 
densate. In our approach, as extensively discussed in ref. [11] , the gluon condensate 
generates the mid range part of the potential, which other authors parametrize as 
being linear. Our calculations indicate that within the approximation 16b , i. e. , 
neglecting higher order condensates, a small reduction in (ft 2 due to the temperature is 
enough to make the existence of ift' impossible. This might be true but might be just 
a consequence of the approximation, which breaks down for very large states. 

There is some uncertainty in the quark masses m c and mj,. One might think that 
a different value of m c would change our results not only quantitatively but also qual- 
itatively. We therefore present in figure 2 plots of the charmonium mass (scaled by 
its mass at zero temperature) as a function of temperature for three different charm 
quark masses m x = 1200, m 2 = 1400 and m 3 = 1640 MeV. For each value of m c the 
octet-singlet energy splitting constant C has to be chosen so as to properly reproduce 
the 1S-2S charmonium mass splitting at T = 0. The chosen values were C\ = —995 
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, C 2 = —885 and C 3 = —737 MeV. We concentrate on the charmonium because it is 
larger than the bottomium and therefore more sensitive to variations of the parame- 
ters. In figure 2a and 2b we show the behaviour with temperature of the charmonium 
fundamental and first excited state respectively. As it can be seen the curves corre- 
sponding to the three quark masses are similar and exhibit the same increasing trend. 
In the same way as charm states have a stronger dependence with temperature than 
bottom states, we observe that lighter charm masses lead to a stronger dependence of 
the bound state with temperature. We have checked that changes in the values of the 
parameters A , b and <p 2 lead to curves with the same aspect of those presented in 
figures 1 and 2. This is also true for the bottomium. From what was said above we 
can see that our results are stable under variations of the parameters. 
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Figure 1: Mass (in vacuum mass units) plotted against temperature (in critical tem- 
perature units). m c = 1640 MeV 
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Figure 2: a) Mass (in vacuum mass units) of the charmonium groundstate as a function 
of temperature (in critical temperature units) for charm quark masses equal to 1200, 
1400 and 1640 MeV. b) The same as a) for the charmonium first excited state. 
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